Experimental implementation of an adiabatic quantum optimization algorithm 
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We report the reahzation of a nuclear magnetic resonance computer with three quantum bits 
that simulates an adiabatic quantum optimization algorithm. Adiabatic quantum algorithms offer 
new insight into how quantum resources can be used to solve hard problems. This experiment 
uses a particularly well suited three quantum bit molecule and was made possible by introducing a 
technique that encodes general instances of the given optimization problem into an easily applicable 
Hamiltonian. Our results indicate an optimal run time of the adiabatic algorithm that agrees well 
with the prediction of a simple decoherence model. 



Since the discovery of Shor'silj] and Grover's al- 
gorithms, the quest of finding new quantum algorithms 
proved a formidable challenge. Recently however, a novel 
algorithm was proposed, using adiabatic evolution 
Despite the uncertainty in its scaling behavior, this al- 
gorithm remains a remarkable discovery because it offers 
new insights into the potential usefulness of quantum re- 
sources for computational tasks. 

Experimental realizations of quantum algorithms in 
the past demonstrated Grover's search algorithm (EMg, 
the Deutsch- Jozsa algorithm 0, S 0| j order-finding- |inl| , 
and Shor's algorithm Recently, Hogg's algorithm 

was implemented using only one computational step ) 
however a demonstration of an adiabatic quantum algo- 
rithm thus far has remained beyond reach. 

Here, we provide the first experimental implementa- 
tion of an adiabatic quantum optimization algorithm us- 
ing three qubits and nuclear magnetic resonance (NMR) 
techniques 0, Q| . NMR techniques are especially at- 
tractive because several tens of qubits may be accessi- 
ble, which is precisely the range that could be crucial in 
determinin g th e scaling behavior of adiabatic quantum 
algorithms [15| . Cornpared to earlier implementations of 
search problems [g, this experiment is a full imple- 
mentation of a true optimization problem, which does 
not require a black box function or ancilla bits. 

This experiment was made possible by overcoming two 
experimental challenges. First, an adiabatic evolution re- 
quires a smoothly varying Hamiltonian over time, but the 
terms of the available Hamiltonian in our system cannot 
be smoothly varied and may even have fixed values. We 
developed a method to approximately smoothly vary a 
Hamiltonian despite the given restrictions by extending 
NMR average Hamiltonian techniques 0|. Second, gen- 
eral instances of the optimization algorithm may require 
the application of Hamiltonians that are not easily acces- 
sible. We developed methods to implement general in- 
stances of a well known classical NP-complete optimiza- 
tion problem given a fixed natural system Hamiltonian. 

We provide a concrete procedure detailing these meth- 



ods. We then apply the results to our optimization prob- 
lem which is known as Maximum Cut or MAXCUT [l^ . 
Our experimental results indicate there exists an optimal 
total running time which can be predicted using a deco- 
herence model that is based on independent stochastic 
relaxation of the spins. 

The evolution of a quantum state during an adiabatic 
quantum algorithm is determined by a slowly varying, 
time-dependent Hamiltonian. Suppose we are given some 
time-dependent Hamiltonian H{t) where < t < T, 
and at t = we start in the ground state of H{0). By 
varying H{t) slowly, the quantum system remains in the 
ground state of H{t) for all < t < T provided the low- 
est two energy eigenvalues of H{t) are never degenerate 
[l8| . Now suppose we can encode an optimization prob- 
lem into H{T). Then the state of the quantum system 
at time t = T represents the solution to the optimization 
problem ^ . The total run-time T of the adiabatic algo- 
rithm scales as where gmm is the minimum separa- 
tion between the lowest two energy eigenvalues of H{t) 
0, El ■ It is the scaling behavior of gmin that will ul- 
timately determine the success of adiabatic quantum al- 
gorithms. Classical simulations of this scaling behavior 
are hard due to the exponentially growing size of Hilbert 
space. In contrast, sufficiently large quantum computers 
could simulate this behavior efficiently. 

Smoothly varying some time-dependent Hamiltonian 
appears straightforward but it contrasts with the tra- 
ditional picture of discrete unitary operations including 
fault tolerant quantum circuit constructions . Fortu- 
nately, we can approximate a smoothly varyin g H amil- 
tonian using methods of quantum simulations j2l| and 
recast adiabatic evolution in terms of unitary operations. 

Discretizing a continuous Hamiltonian is a straightfor- 
ward process and changes the run time T of the adia- 
batic algorithm only polynomially |l9l |. For simplicity, 
let the discrete time Hamiltonian H[m\ be a linear inter- 
polation from some beginning Hamiltonian H[0\ — Hb to 
some final problem Hamiltonian H[M] — Hp such that 
H[m] = {m/M)Hp 4- (1 - m/M)Hb. The unitary evolu- 
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tion of the discrete algorithm can be written as: 



-i{{l-m/M)Hb + {m/M)Hp)At 



(1) 



where At = T/{M + 1), and M + 1 is the total number 
of discretization steps. The adiabatic limit is achieved 
when both T, M oo and At 0. 

Full control over the strength of Hi, and Hp is needed 
to implement Eq.(^. However, this may not necessarily 
be a realistic experimental assumption. We will next 
show how the discrete time adiabatic algorithm can still 
be implemented when Ht and Hp cannot both be applied 
simultaneously and when they are both fixed in strength. 

When both Hi, and Hp are fixed, we can approxi- 
mate Um to second order by using the Trotter formula 
exp((^ + B)At) = ex.p{AAt/2)exp{BAt)exp{AAt/2) + 
©(At^) [23. Higher order approximations can be con- 
structed if more accuracy is required. 

Now suppose Hh and Hp are both constant. Since any 
unitary matrix is generated by an action —iHAt, we can 
increase the effect of a constant Hamiltonian H by length- 
ening the time At. Thus, we can implicitly increase the 
strength of Hi, and Hp even when they are constant by 
simply increasing the time during which they are applied. 

This technique also allows cases when the accessible 
Hamiltonians are not of the required strength, for exam- 
ple when we are given iJ^ = gHi, and Hp — hHp but still 
wish to implement Hj, and Hp . Using all of the described 
techniques, we can now write Um as: 



U„ 



'iH'Ml-,n/M)At/2g] 



-iH'[{m/M)At/h] 



(2) 



where AoB = ABA. Each discretization step is of length 
(1 - m/M)At/g + {■m/M)At/h, which is not constant 
when g ^ h. As an illustration consider Fig. 

In this experiment we choose At = T/{M +1) to be 
constant as we vary the number of discretization steps 
M + 1. This way, the total run time T increases with 
M + 1, allowing us to test the behavior of the algorithm 
when approaching one of the conditions for the adiabatic 
limit. Even when the discrete approximation is not close 
to the adiabatic limit, the implemented algorithm can of- 
ten find solutions using relatively few steps but lacks the 
guaranteed performance of the adiabatic theorem ji^l . 

Adiabatic evolution has been proposed to solve gen- 
eral optimization problems, including NP-complete ones. 
In this general setting, the algorithm can depend on the 
existence of a black box function or the usage of large 
amounts of workspace. Our goal here is to optimize a 
hard natural problem in a way that avoids these difficul- 
ties. We will first describe which problem we chose and 
later on explain why it does not require ancilla qubits. 

We found the MAXCUT problem to be a well-suited 
problem to demonstrate an adiabatic quantum algorithm 
because it allows a variety of interesting test cases. It 
also has applications in the study of spin glasses and 




FIG. 1; (a) Illustration of Eq.l(2J. The shaded and clear boxes 
denote the strength and duration of the Hamiltonians Hb and 
Hp respectively, (b) Illustration of a graph consisting of 3 
nodes and 3 edges. The edges carry weights W12 , W13 , and W23 . 
When min(wi-, ) — W23 as indicated by the length of the edges, 
the MAXCUT corresponds to the drawn cut. The solution is 
therefore s = 100 and also s = Oil due to symmetry. This 
symmetry can be broken by assigning the weights wi, W2, and 
W3 to the nodes. 



VLSI design [2J|, among others. The decision variant 
of the MAXCUT problem is part of the core NP-complete 
problems , and even the approximation within a fac- 
tor of 1.0624 of the perfect solution is NP-complete '2^. 

The MAXCUT problem can be understood as follows. 
A cut is defined as the partitioning of an undirected n- 
node graph with edge weights into two sets. We define 
the payoff as the sum of weights of edges crossing the 
cut. The maximum cut is a cut that maximizes this 
payoff. By assigning either Si — or = 1 to each 
node i, depending on its location with respect to the cut, 
the MAXCUT problem can be restated as finding the n 
bit number s that maximizes the payoff. An extension of 
the MAXCUT problem is to let the nodes themselves carry 
weights, which can be regarded as the nodes having a 
preference on their location. As an illustration consider 
a graph with three nodes as drawn in Fig. 

The payoff as a function of the cut defined by s is given 

by 



Pis) 



i,3 



.(1 



(3) 



where Wij are the edge weights, Wi denotes the preference 
of the nodes to be on the 1 side of the cut, and is the 
value of the i-lh bit of s, for < s < 2" — 1. 

The smallest meaningful test case of the MAXCUT prob- 
lem requires 3 nodes and admits a variety of interesting 
cases by varying Wi and Wij . We aimed at two goals when 
choosing a representative set of weights. First, we wanted 
the minimum energy gap gmin to be smaller than the one 
for a 3-qubit adiabatic Grover search. Second, we wanted 
a resulting energy landscape with both a global and lo- 
cal maximum such that a greedy classical search would 
incorrectly find the local maximum half the time [2^ . 
These goals are met by the choice wi = W2 = = 2, 
W12 = 2, W13 = 1, W23 = 3. The payoff function for 
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this set of weights is P(s) = [0 677598 6] where s = 
[000 001 010 Oil 100 101 110 111]. The global maximum 
lies at s = 101 so the answer on the quantum computer 
following measurement should be |101), and not at the 
local maximum s = 110. 

In the quantum setting, this payoff function P{s) can 
be encoded into the Hamiltonian Hp by rewriting Eq. 10) 
using Pauli matrices: 

Hp = J2MI -'^zi)/2 + Y,W'^jiI -'^^^<^z3)/'2 (4) 
i i<j 

where / is the 2"x2" identity matrix and azi is the Pauli 
Z matrix on spin i. The identity matrices in the equa- 
tion above only lead to an overall phase which cannot be 
observed, and hence they can be ignored. The diagonal 
values of Eq.Q are equal to P(s). Because of the direct 
encoding of -P(s) into Hp no black box function or ancilla 
qubits are required, which makes this a full implementa- 
tion of an optimization problem. 

Similar to Eq. I^J , the natural Hamiltonian of n weakly 
coupled spin-1/2 nuclei subject to a static magnetic field 
Bq is well approximated by 

H = — ^01^0-21/2 -I- ^ TrJijaziCrzj/2 + Henv (5) 
i i<j 

where the first term represents the Larmor precession of 
each spin i about —Bq, and is its Larmor frequency. 
The second term describes the scalar spin-spin coupling 
of strength between spins i and j. The last term 
represents coupling to the environment, causing decoher- 
ence. Note the resemblances between Tl and Hp. 

Despite the similarities, the spin-spin couplings of 
Eq. (|3J) are generally different from a randomly chosen set 
of weights. Therefore, we require a procedure to turn the 
fixed Jij into any specified weights Wij . This is achieved 
using refocusing schemes that are typically used to turn 
on only one of the couplings while turning all others off 
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We have modified a refocusing scheme to effectively 
change the couplings to any arbitrary value. Consider the 
pulse sequence drawn in Fig. |3 Based on this scheme, 
we can derive the under-constrained system {a + f3 — j — 
S)Ji2 = ^12, (a - /3 - 7 -I- (5) Ji3 = Wis, and {a - (3 + 
7 — (5) J23 = 1^23, which can be solved for positive a, /3, 
7, and d such that Jy —>■ Wij. 

The single weights Wi are implemented by introducing 
a reference frame for each spin i which rotates about —Bq 
at frequency {tOi ~ Wi)/2. In order to apply the single 
qubit rotations of our refocusing scheme on resonance, 
we apply the reference frequency shift only during the 
delay segment a, which we can always choose to be a 
positive value. Thus, Hp is implemented by applying the 
refocusing scheme from Fig. |21 while going off-resonance 
during the delay segment a. 

A full implementation of an adiabatic algorithm also 
requires a proper choice of Ht. We choose H), = a^i 




FIG. 2: Refocusing scheme to effectively change Jij into Wij. 
The horizontal lines denote qubits 1, 2, and 3 and time goes 
from left to right. The black rectangles represent 180° rota- 
tions. The delay segments are of length q, /9, 7, and S. When 
all segments are of equal length, all couplings are effectively 
turned off because axie~^''^^'^'^^axi = e"'^'°'^3'. In our 
experiment, a = 0.42 ms, /3 — ms, 7 = 4 ms, and S — 2.9 
ms in the last slice M + 1. The RF pulses that implement 
Hb' perform 33.75° rotations on the qubits in the first slice. 



for several reasons. First, its highest two excited states 
are non-degenerate. Second, it can be easily generated 
using single qubit rotations, and third, its highest excited 
state is created from a pure state with all qubits in the 
jo) state by applying a Hadamard gate on all qubits (we 
require the initial state to be the highest excited state of 
Hh because we are optimizing for the maximum value of 
Hp). 

The full adiabatic quantum algorithm is now imple- 
mented by first creating the highest excited state of Hi,. 
We then apply M+1 unitary matrices as given by Eq.Q 
and illustrated by Fig^. Accordingly, from slice to slice, 
we decrease the time during which Hj, is active while in- 
creasing the time during which Hp is active. Finally, we 
measure the quantum system and read-out the answer. 

We selected ^'^C-labeled CHFBr2 for our experiments. 
The Hamiltonian of the ^H-^^F-^^C system is of the 
form of Eq.llsj with measured couplings Jhc = 224 Hz, 
Jhf — 50 Hz, and Jpc — —311 Hz. The interaction 
with the Br nuclei is averaged out, contributing only to 
Tienv ■ Experiments were carried out at MIT using an 11.7 
Tesla Oxford Instruments magnet and a Varian Unity In- 
ova spectrometer with a triple resonance (H-F-X) probe 
from Nalorac. 

The experiments were performed at room temperature 
at which the thermal equilibrium state is highly mixed 
and cannot be turned into the required initial state by 
just unitary transforms. We thus first created an approx- 
imate effective pure state as in ref. |^ by summing over 
three temporal labeling experiments. 

In our experiments, we actually implemented O.bHp 
and 0.5887Hi, instead of Hp and Hi,. This ensures that 
the error due to the 2nd order Trotter approximation is 
bounded by je^j^ < 0.0356 where is the contribu- 

tion of the i-th undesired Pauli matrix. We also choose 
g so the applied RF field does not heat the sample, and 
g ^ h so Jij can be ignored when applying Hi,. All of 
these choices result in a total experimental time that is 
within the shortest T2 decoherence time 6]. We recon- 
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structed the traceless deviation density matrices upon 
completion of the experiments using quantum state to- 
mography 6J. 

We executed this algorithm for several M (with Wi 
and Wij as listed above Eq.Q). Since we chose At to 
be constant, this meant increasing the run-time T of the 
algorithm. The reconstructed deviation density matrices 
are shown in Fig. 13 The plots clearly display the ex- 
pected pure state |101). The local maximum at s — 110 
has a decreasingly small probability of being measured 
for increasing M. Simulations using Eq.© show that 
this optimization algorithm performs better for increas- 
ing M . We wanted to verify whether this is indeed true 
experimentally. 
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FIG. 4: Experimental performance of the adiabatic algorithm, 
(a) Plot of the error as a function of M. The error measure is 
the trace distance D{p, o) = |p — a|/2 where a is the traceless 
deviation density matrix for M = 400, approximating M —f 
oo, and p equals the ideal expected (o), the experimentally 
obtained ( x ) , or the ideal expected traceless deviation density 
matrix with decoherence effects (o) W^. The minimum error 
occurs at about M = 60 indicating an optimal run-time of 
the algorithm, (b) A similar observation can be made when 
plotting 1 101} (101 1 as a function of M. 
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FIG. 3: Plot of the absolute value of the deviation density 
matrix for M = 100 (T = 374 ms), M = 30 (T = 115 ms), 
and M = 15 {T — 59.2 ms), adjusted by an identity portion 
such that the minimum diagonal value equals zero. The scale 
is arbitrary but the same for each plot. 

For this purpose, we estimate the error of our obtained 
deviation density matrices compared with the ideal case 
of M — oo. Fig. plots the trace distance as a func- 
tion of M, using the same arbitrary scale as in Fig. O 
From the plot, we observe there exists an optimal run- 
time of the algorithm, corresponding to 0.226 seconds 
in our experiment. This optimal run time is in good 
agreement with the prediction of a previously developed 
simple decoherence model . Predicting the impact of 
decoherence has already provided invaluable insight into 
estimating errors in previous experiments and we 
believe continued effort towards understanding decoher- 
ence will greatly benefit experimental investigations of 
quantum systems. 

In conclusion, we have provided the first experimen- 
tal demonstration of an adiabatic quantum optimization 
algorithm. We show a concrete procedure turning a con- 
tinuous time adiabatic quantum algorithm into a discrete 
time version, even when certain restrictions apply to the 
accessible Hamiltonians. Our results indicate that there 
exists an optimal run-time of the algorithm which can be 
roughly predicted using a simple decoherence model. We 
believe this implementation opens the door to a variety 
of interesting experimental demonstrations and investi- 
gations of adiabatic quantum algorithms. 
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